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Electron collisions for a two dimensional Fermi liquid (FL) are shown to give a quasiparticle 
damping with interesting frequency and temperature variations in the BCS superconducting state. 
The spin susceptibility which determines the structure of the damping is analyzed in the normal 
state for a Hubbard model with a constant on-site Coulomb repulsion. This is then generalized 
to the superconducting state by including coherence factors and self energy and vertex corrections. 
Calculations of the NMR relaxation rate reveal that the FL damping structure can reduce the 
• Hebel-Slichter peak, in agreement with data on the organic superconductor (MDT-TTF)2Aul2. 

I However, the strongly suppressed FL damping in the superconducting state does not eliminate the 

Hebel-Slichter peak, and thus suggests that other mechanisms are needed to explain the NMR data 
on (TMTSF)2C104, the BEDT organic compounds, and cuprate superconductors. Predictions of 
Qj^, the temperature variation of the damping and the spin response are given over a wide frequency 

■ range as a guide to experimental probes of the symmetry of the superconducting pairs. 
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0^ ' The damping of quasiparticles via electron collisions is expected to produce a variation of the resistivity in a 
' standard Fermi liquid described by a spherical Fermi surface. In ordinary metals this contribution is so weak that it 
is hardly detectable 0. However, recent discoveries of a dominant resistivity contribution in several anisotropic 
metals such as organic superconductors and layered cuprate and sulfide alloys provide an interesting challenge from 
the theoretical point of view since the observed magnitudes of such resistivities are very large. Examples of the 
unconventional resistivities are shown in Figure |l[ 
^ , In view of the conventional wisdom regarding the expected weakness of the electron-electron scattering contribution 
'"^ to the resistivity, alternate interpretations of the TiS2 data were proposed on the basis of unusual phonon scattering. 
^ However, decisive evidence for the Fermi liquid origin of the scattering was discovered by Julien jsj in the infrared 
O spectra of TiS2. These spectra show a remarkable frequency and temperature variation in accordance with the 
Luttinger damping for a Fermi sphere that was derived to all orders of perturbation theory for the Coulomb interaction 
I> [Q. The huge value of the resistivity of TiS2 is on par with the data for the Nd2-xCexCu04 cuprate which also has a 
' layered structure. It is also comparable to the resistivity of many organic metals such as that of (TMTSF)2PF6 shown 
, in Figure ||. The organics have quasi- two dimensional Fermi surfaces with a high degree of anisotropy. By contrast, 
■ the much smaller resistivity of lead shown in Figure ^ provides an example of the scale set by strong electron-phonon 
contributions. 

The purpose of the present work is to calculate the quasiparticle damping for a two dimensional electron system due 
to spin fluctuation scattering arising from the constant on-site Coulomb interaction U in a Hubbard model. We first 
compute the standard normal state damping and then focus on the influence of an isotropic superconducting energy 
gap on the frequency and temperature variation of the self energy. The superconducting energy gap is expected to 
reduce the damping in a characteristic way that depends on the symmetry of the order parameter as well as the 
source of the damping. The spin susceptibility, which enters in the calculation of the damping, is first analyzed for 
the normal electrons and then generalized to include the BCS coherence factors in the superconducting state. Our 
goal is to use the computed damping in a calculation of the NMR relaxation rate, with particular interest in the 
superconducting state response. The relaxation rate is obtained from a momentum average of the spin susceptibility 
that includes the calculated FL self-energy corrections. Vertex corrections to the susceptibility are also calculated 
and shown to be small in the regime considered in the NMR studies. 

We apply our results to experimental NMR measurements on the organic superconductors. These provide insight 
into the structure caused by the damping as well as the symmetry of the superconducting state. The modification of 
the Hebel-Slichter (HS) peak is found to be particularly instructive in the case of the anomalies seen in the organics, 
although the general phenomenon is relevant to the high temperature cuprate superconductors as well. 



1 



Luttinger |j] originally derived the classic three dimensional Fermi liquid property of a quasiparticle damping arising 
from electron collisions that vanishes at the Fermi energy at zero temperature and then follows a quadratic variation 
in temperature T and frequency uj (measured from the Fermi energy) . Hodges et al. Q and others Q calculated the 
quasiparticle damping and resistivity for a cylindrical Fermi surface, which introduces logarithmic corrections to the 
damping at low T and w. 

We first calculate the two dimensional Fermi liquid damping in the normal state by performing a momentum average 
over the exact spin susceptibility for a non-interacting electron gas. We then proceed to compute the susceptibility as 
well as the self energy in the superconducting state, and probe the influence of the calculated damping on the NMR 
response which involves a momentum-averaged spin susceptibility. 

The original motivation for our present work was the mysterious absence of a Hebel-Slichter peak in the NMR 
response of high temperature cuprate superconductors |^]. The resistivity of the optimally doped cuprates above 
the superconducting transition is typically linear in temperature in contrast to the Fermi liquid behavior considered 
here. A similar physical origin of electron-electron scattering in both cases may be plausible, because alloying changes 
the linear T to a variation in many of the cuprates |8|. Historically the HS peak in the NMR relaxation is 
considered to be one of the key successes of the BCS theory since it occurs in various ordinary metals The 
organic superconductors yield cases where the HS peak appears, (MDT-TTF)2Aul2 jlOll , and also examples like 
(TMTSF)2C104 [jlll where the peak is suppressed. The NMR relaxation in the BEDT [[12[| compounds resembles the 
cuprate anomalies, such as a behavior at low T in the superconducting state. Hence the present analysis is directed 
at the organics which often have the quadratic T variation of the resistivity that is compatible with our Fermi liquid 
analysis. 

The absence of the coherence peak has attracted considerable previous theoretical interest. Hasegawa and Fukuyama 
p3| derived the NMR response for s-wave and other symmetries of the energy gap appropriate to organic supercon- 
ductors. They obtained weaker, but nonetheless finite, HS peak structure for both singlet and triplet pairing states in 
which the order parameter exhibits line nodes on the Fermi surface. Their calculation did not include the quasiparticle 
damping. 

Other groups have examined the effects of a strong enhancement of the quasiparticle damping near Tc, which 
broadens the BCS singularity in the superconducting density of states. Some theory groups have invoked a variety 
of phcnomenological models for the damping as a function of temperature . Typically, these authors use a power 
law variation for the T dependence of the damping (with a behavior as one example) and neglect the frequency 
variation. Phcnomenological models have been proposed to fit the observed spin dynamics in the cuprates ||15| . The 
"marginal" Fermi liquid response hypothesis has been used to generate a damping model which has also been applied 
to the NMR measurements |l6|| . 

Tight binding energy band models have also been investigated in the context of spin dynamics of cuprates by Bulut 
and Scalapino |17| and by Levin's group ||l8|. The NMR relaxation for s and d-wave energy gaps has been calculated 
p^ , pT| using phenomenological models for the damping. The quasiparticle damping from spin fluctuations in the 
superconducting state has also been independently computed within a tight binding model [T9| ] . 

Our analysis of the influence of the damping on the NMR spectra in the superconducting state is based on a standard 
approximation for the momentum-averaged spin susceptibility. This allows analytic expressions to be derived and 
used in the calculation of the damping, thus reducing the number of numerical integrations. However, this limits the 
present approach to reasonably isotropic Fermi surfaces, and cannot account for features such as nesting. 

Phonon damping has furnished yet another more traditional explanation for the suppressed HS peak, provided 
that a large electron-phonon coupling combined with a small Coulomb pseudopotential are assumed Phonon 
contributions to the resistivity may be one distinguishing feature of the mechanism responsible for the NMR anomalies, 
and the organic resistivity data in Figure |l| provide a challenge in this regard since the variation is difficult to 
reconcile with phonon scattering. 

Considering the wide array of theoretical proposals for the NMR anomalies, it is important to make distinctions 
based on specific features such as the response over a wide range of temperature, and consistency of parameters with 
other experimental clues. The damping which characterizes the organic superconductors considered here is one 
example. The temperature dependence of the NMR relaxation rate at low T is also a primary constraint on theory 
since it is particularly sensitive to the symmetry of the order parameter, while the HS peak structure is sensitive to the 
damping as well as the gap symmetry. Thus our calculations provide insight into the role of damping in the organic 
superconductors and the NCCO cuprate which also exhibits a resistivity above the superconducting transition. 

A theoretical basis for the anomalous linear T variation of the resistivity that is an ubiquitous feature in all of 
the optimally doped high temperature superconductors discovered so far is the nested Fermi Liquid theory (NFL) 
|21|. The corresponding calculations for the NFL damping in the presence of an isotropic energy gap by Rieck et al. 
1 22 1 reveal a very dramatic suppression of the nested spin susceptibility in the superconducting state which strongly 
reduces available scattering states for electron collisions. Hence the NFL damping is greatly reduced at frequencies 
lower than thrice the energy gap A. The latter damping structure is compatible with microwave surface resistance 
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data on the YBCO superconductor in the vicinity of Tc where the damping drops by four orders of magnitude within 
a few degrees of T^- There are also numerous earlier calculations that found a sharp suppression of the damping 
in the case of phenomenological models that are extensions of the "marginal" Fermi liquid hypothesis for the 
susceptibility. These are discussed in reference [|2|. 

Recently Won and Maki have shown that a nesting model yields NMR relaxation without a HS peak if indeed 
the scattering processes near the nesting vector form the dominant contribution to the susceptibility. 

We develop the formalism for the susceptibility and self energy in section II, and present the calculated Fermi liquid 
damping results. The NMR relaxation is discussed in section III, and the conclusions of our study are in section IV. 



II. FORMALISM 

We consider the Hubbard Hamiltonian 

H^Yl ^k4,,Ck,, + ^ H 4+q,TS,T4-qaCk,i> (1) 
k,cr p,q,k 

where ek is the energy of an electron in two dimensions, and cj^ ^ and Cj^ ^ are the electron creation and destruction 
operators. The constant on-site Coulomb repulsion U restricts the scattering to spin fluctuations. Within the Born 
approximation considered here, the dominant effect is from the non-interacting spin susceptibility and the cross section 
is proportional to . In other words the self energy requires a momentum sum over C/^xIQ:"-^)- Higher order RPA 
corrections will naturally enhance the cross section, but are not included here. 



A. Normal State 



The spin susceptibility in the normal state for non-interacting electrons is 

/ (ek+q) - / (fk) 



Lo - ek^ 



Ek + »r 



(2) 



where /(e) is the Fermi function. Within the leading order Born approximation, the imaginary part of the self energy 
arising from electron-electron scattering is 
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where x" is the imaginary part of x{q^,Lu). 

We first verified numerically that the imaginary part of the susceptibility in Eq. (H) does not vary significantly with 
temperature for a normal two-dimensional electron gas, and hence one may use the zero temperature analytic result 
in calculating F. Defining the frequency w in units of hkp/m and the wavevector q in units of the Fermi wavevector 
kp, the T = analytic form for x" is 



X" (q,^) = 



N{0) 



q 2 



q 2 



(4) 



where -/V(O) is the density of states at the Fermi level, and 9{x) is the Heaviside function that is unity for a; > and 
zero otherwise. Hodges ^ and others |^ have used the above analytic form to calculate the damping in the asymptotic 
limits of small T and to, by performing the momentum integrations (Eq. ^) to obtain the transport lifetime. Note 
that the imaginary part of the self energy is F = [1 — J{uj)\^^'h/T where t is the quasiparticle transport lifetime [ p5| . 

To obtain a damping over the entire range of frequency and temperature, we employ an alternate approximation: 
the susceptibility in Eq. (|^) is first averaged over momentum and gives at low frequencies (cj < ep), 



(X"(q,c.)>. 



7rA^(0) 
W 



-UJ 



(5) 



where W is the energy bandwidth, and iV(0) is the density of states at the Fermi energy. We note that the momentum- 
averaged susceptibility is independent of temperature for the Fermi liquid. Replacing this average in Eq. (0), we 
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evaluate T numerically and find a quasiparticle damping that is quadratic in frequency and temperature. The limiting 
cases are 

r{uj = 0) = BT^ (6) 

r(r = 0) = Clj^ (7) 

where B = U'^N'^{Q)ttIJ2W and C = U"^ N'^ {Q)tt / 2W . The present method does not yield the T^lnT dependence 
found asymptotically [§j6|, but the InT correction is almost indistinguishable from at low T. 

The quadratic temperature variation of our computed FL damping at two different frequencies is shown in Figure |^ 
for a coupling UN{Q) — 1 and a bandwidth W — 1 eV. These parameters are used throughout this work. As expected, 
the value of the FL damping seen in Figure ^ is quite small. For comparison, the experiments that we examine suggest 
much larger estimates of the damping which may perhaps indicate anisotropic Fermi surface contributions, RPA 
enhancement of the spin susceptibility, or higher order scattering. Also, the Hubbard model restricts the Coulomb 
interaction to a point in real space and thus samples only the zero angular momentum scattering channel. 



B. Superconducting State 



In the absence of the self energy corrections from particle collisions, the BCS susceptibility in the superconducting 
state is modified by the presence of coherence factors due to an energy gap, and thus becomes 
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k,q- 
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where the coherence factors are 
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^ _|_ ekCk+q + AkAk+q 
^^k-Ek+q 



(8) 
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A| the 



and Ak is the superconducting energy gap, Ck the energy dispersion in the normal state, and E^ = 
quasiparticle energy in the superconducting state. All the calculations described in this article are carried out using 
a superconducting gap of s-wave symmetry. The temperature dependence of the gap is taken to have the standard 
form that fits the solution of the BCS weak coupling gap equation. 



A(r) = Aotanh 1^1.76-^^ -ij 

where Aq = l.T&kBTc- In the limit of negligible damping, (5^0+, the above susceptibility (Eq. 
sharp Hebel-Slichter peak in the NMR response. 

The quasiparticle damping due to spin fluctuations in the superconducting state is pql 
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gives the standard 
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where n{Lo) and f{Lu) are the Bose and Fermi functions respectively. To simplify the multiple integrations needed 
for the damping, we note that the momentum variation of the susceptibility arising from a free electron dispersion 
is relatively smooth, in contrast to tight binding models considered, for example, by Quinlan et al. fl^ . Hence we 
approximate the damping calculation by taking a momentum average of the susceptibility (for uj > Q) 
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The resulting calculated susceptibility average is shown in Figure |3|, where the peak just below Tc at low frequencies 
(solid line) is a consequence of the coherence factors. At higher frequencies, as seen in the case of w = O.SAq (dashed 
line), the susceptibility drops off rapidly as T decreases though there is no coherence peak below Tc- This characteristic 
drop of the susceptibility below Tc is the key to determining the temperature variation of the quasiparticle damping 
in the superconducting state. We have used here a superconducting transition temperature Tc — 4.2 K, which 
corresponds to the organic compound (MDT-TTF)2Aul2. 

We next compute the quasiparticle damping in Eq. ( pT[ ) as a function of frequency and temperature for a super- 
conductor with Tc = 4.2 K. The frequency dependence of F is shown in Figure ^ at three different temperatures, and 
it reveals the dramatic drop in F at low lu caused by the isotropic energy gap. Note that at zero temperature there is 
no structure in F below SAq. At higher temperatures, F displays an unusual frequency variation below 3A(r), and 
displays a roughly quadratic increase at higher frequencies. 

The temperature dependence of the computed F is shown in Figure |^ at three different frequencies. For frequencies 
below Lj = 3Ao, the damping drops to zero at low T, but reduces to a finite value even in the zero temperature limit 
when the frequency exceeds thrice the energy gap. 

For the calculation of the NMR relaxation rate described in the next section, we reduce the computational complexity 
by developing a model fit to the numerically calculated frequency and temperature dependence of the Fermi Liquid 
F at a given value of Tc- Above Tc, the damping is presumed quadratic in frequency and temperature. 



III. NUCLEAR RELAXATION RATE 



The standard expression for the NMR relaxation in normal metals is 

J- = \A\'J-Y.x"iq,^o) (13) 
TiT ujQ ^ 

where |^| denotes the hyperfine coupling and cjo ~ 10 MHz is the typical radio frequency of the measurement. For a 
Fermi Liquid, the momentum average of the spin susceptibility in Eq. (|^) is linear in frequency, and thus produces 
the familiar Korringa relation for the nuclear relaxation that is of the form 

^^|Ap7V2(0) = constant. (14) 

If the self energy is only weakly momentum dependent, the general form of the susceptibility in Eq. (^3|) can be 
averaged over k and q independently, and the Korringa behavior persists with very little contribution from the damping 
in the normal state. However, an isotropic energy gap produces a divergent density of states in the superconducting 
state which is naturally quite sensitive to the form of the damping. Since the above analysis demonstrates that the 
structure of the damping is particularly important at frequencies comparable to the energy gap, these features are 
important for the NMR spectra in the superconducting regime. 



A. Self Energy Corrections 

The NMR relaxation rate is computed using the form |27) , 

^ = 2 1^1^ N\^) r r _^^EL= + f{E) [1 - J{E')] . (15) 

where A'^(O) is the normal state density of states at the Fermi energy, and the spin hyperfine coupling \A\ is taken to 
be constant, thereby neglecting crystalline anisotropy. The delta function in the original expression for x" in Eq- ( |l^ ) 
is replaced by a Lorentzian of width F, the quasiparticle damping computed in Eq. (|l l| ) . Thus the conventional BCS 
result is recovered in the limit F ^ 0. The self energy corrections included in Eq. ( |15| ) depend on both temperature 
and energy, and enter as F = i [F(£',r) + T{E',T)]. In all the calculations discussed below, the damping is defined 
as 

F = aFpL (16) 

where a is varied to show the consequences of enhancements beyond the Born approximation used here with one 
standard set of parameter values W = 1 eV and UN{0) = 1. A remarkable feature of the organic compound 



5 



resistivities shown in Figure |^ is that their high values would indicate very short mean free paths if a standard 
transport model is applied. However, the temperature variation of the resistance is compatible with Fermi liquid 
behavior, so that the enhancement remains anomalous in these cases even though the likely physical origin is electron- 
electron scattering. Anisotropic Fermi surface effects or higher order scattering such as the RPA corrections may be 
likely suspects for this mystery. 

The relaxation rate data of Takahashi et al. [|l0| on the organic superconductor (MDT-TTF)2Aul2 is shown in Figure 
^ along with the theoretical curves calculated using Eg . (|l5|). The solid curve is calculated using the superconducting 
state quasiparticle damping computed above (Eq. |ll| ) for the Fermi liquid. To examine the influence of the energy 
gap on the relaxation rate, we compare this to the dot-dashed curve calculated using the normal-state damping that 
is quadratic in frequency and temperature at all T and uj. Note that in both these cases the Hebel Slichter peak is 
strongly reduced over the BCS result (dashed line) obtained in the absence of a broadening F. However, in order 
to fit the (MDT-TTF)2Aul2 data, the amplitude enhancement of the calculated Fermi Liquid damping needs to be 
a = 40. The standard case with a = 1 resembles the dashed BCS curve because the damping is very small. 

The case of (TMTSF)2C104 presents an interesting challenge since the NMR data of Takigawa et al. clearly 
show the absence of a Hebel-Slichter peak as seen in Figure n. The resistivity of this organic superconductor has a 
temperature dependence that is close to in some samples p8| , while carefully quenched samples exhibit a linear T 
resistivity |29|. This material is notable for nesting of the Fermi surface which gives rise to spin density wave (SDW) 
phase transitions in a magnetic field as discussed in terms of the nested orbit quantization by Gorkov and Lebed . 
Our goal is to fit the NMR data using the Fermi liquid damping to see whether the damping alone can account for the 
suppression of the Hebel-Slichter peak. The results shown in Figure ^ reveal that the Fermi liquid damping actually 
gives a weak, but nevertheless significant peak in the NMR relaxation just below even for very large damping cross 
sections. The two curves in Fig. ^ are calculated using an increase in amplitude of the damping by factors of a = 
40 and 400 respectively, and both curves are normalized to the value of 1/TiT at Tc = 1.06 K. The a = 400 curve 
resembles the data, indicating that an anomalously large damping may indeed be the cause of the suppression of the 
HS peak. 

The data shown in Figure || provide another illustration of non-Fermi Liquid behavior in the organic compound 
K-(BEDT-TTF)2CuN(CN)2Br. The NMR relaxation rate data of Kanoda et al. ^ on this organic compound are 
shown (dots) along with the theoretical curve (solid line) calculated using the Fermi Liquid damping and enhancement 
a = 40. Unlike the calculated curve, the data show no evidence of a Hebel-Slichter peak, and also drop much 
more sharply at low temperatures than the calculated curves. The conventional BCS result with a more pronounced 
coherence peak is also shown (dashed line) . The data for this organic compound exhibit a temperature dependence 
at low T, which is suggestive of an unconventional pairing state with line nodes on the Fermi surface, e.g. d-wave 
pairing. In the normal state the data are also anomalous in that they do not follow a Korringa law | ]3l[ |. A peak in 
1/Ti is observed at 50 K which tends to vanish under pressure This pressure variation of TiT in the normal 

state may arise from spin fluctuations whose contributions increase with the degree of nesting of the Fermi surface. 



B. Vertex Corrections 



The NMR relaxation rate discussed above has been calculated in the presence of self energy corrections to the 
spin susceptibility. Schrieffer has raised the issue of the importance of vertex corrections for pairing schemes 
that rely on the exchange of spin fluctuations, suggesting that the pairing coupling may be strongly suppressed. We 
have addressed this question in the case of the lowest order bubble susceptibility ^(q', w). We are interested in vertex 
corrections to the NMR relaxation rate consistent with the Ward identity. Since 1/TiT = limi^^o x"(9i'^)/'^i it 
suffices to calculate the correction to the momentum averaged imaginary part of the spin susceptibility. Wermbter and 
Tewordt have used strong coupling Eliashberg theory to compute the vertex corrections within a 2D Hubbard 
model for an interaction arising from the exchange of spin and charge fiuctuations as well as phonons. In their 
calculations, they approximate the bubble susceptibility xiQj by its average over momentum. Their approach is 
thus a good starting point for our calculation which is carried out in the weak coupling limit for the lowest order 
electron- hole bubble. 

In the absence of self-energy and vertex corrections, the momentum average of the imaginary part of the suscep- 
tibility is given by the expression for 1/TiT in Eq. (^a). Including vertex corrections to order U^, the one-loop 
contribution relevant to the spin fiuctuations becomes p4[ 

X"{u;) =TT dE ^ ' 



VE^~A^y/{E + u;)^-A^ 
[f{E)[l - Ji{E,Lo) + tt''J2{E,u;)] - f{E + uj)[l - Ji{E,lu) + 7rV^(£;, c^)]] (17) 
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where = {x" {(\, uj)) q is the momentum-averaged susceptibility. The correction terms Ji, J2 are given by 

/(-pi) + n(i-) [~ dn 
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where f{x), n{x) are the Fermi and Bose functions respectively, and 

U^N^{0) . _ „ . \E + uj-iy\\E-iy\+A^ 



J2{E,UJ) 
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(18) 



(19) 



The term in Eq. ( |17| ) is obtained from J2 by substituting ]{v — E) for f{v — E — u). 

The NMR relaxation rate is obtained from the lu hmit of x"(w) in Eq. (p^). We therefore calculate Ji(i?,0) 
and J2(£',0), noting that J'2(E,Q) = J2{E,0). The /i2 integration in Eq. (|l|) can be carried out analytically, giving 
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In the low temperature limit, T — > 0, we can obtain an estimate of the vertex contribution by replacing coth (x/2T), 
and tanh(x/2T) by sgn(x), and the /ii integration in Eq. ( po|) can also be carried out analytically. The remaining v 
integrations in Ji and J2 (Eq. |l^) are carried out numerically. In the limit cj — > 0, Eq. (FtI) becomes 



-i- OC \iin -x"{uj) 

1 il u]—>0 to 
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(21) 



The resulting total vertex correction, —Ji{E,0) + tt'^J2{E,0) is shown in Figure ^. We find that the contribution of 
Ji is approximately an order of magnitude larger than J2 for T < Tc. Since Ji is proportional to A^, the correction is 
larger for a system with higher Tc- This is illustrated in Figure ^ for Tc = 10, 20, 100 K. We find that Ji decreases with 
increasing temperature, and vanishes as T ^ Tc since A ^ 0. In this limit, J2(i?, 0) —U^N^{0)[tt^T^ + E^]/4W^, 
which is small for E <C W. The temperature dependence of J2 reveals that J2{T < Tc) < J2{Tc) for all E. Thus the 
total correction term decreases as T increases. The transition temperatures of the organics are at most of the order 
of 10 K, and the total vertex correction for Tc = 10 K is less than 3%. We also find from Figure ^ that at a fixed 
temperature, regardless of the value of Tc, the correction does not depend strongly on E. It can thus be taken out of 
the integral in Eq. ( pl| ) and the net result is a renormalized value of the interaction strength U. 

We basically concur with the conclusions of Wermbter and Tewordt who estimated vertex contributions by 
another method and found them to be negligible for their choice of parameters. 



IV. CONCLUSIONS 



We have calculated the quasiparticle damping F arising from electron-electron collisions in a two dimensional 
Fermi Liquid. We have presented results for the temperature and frequency variation of F in the normal and the 
superconducting state with an isotropic energy gap. The damping is found to be quadratic in temperature and 
frequency above the superconducting transition, consistent with the T^ resistivity of the organic superconductors and 
the Nd2-xCexCu04 cuprate. 

The computed damping follows an unusual frequency variation in the superconducting state for frequencies 
uj < 3A(T). The reduction in available scattering states due to the opening of an isotropic energy gap causes 
the quasiparticle damping to drop at temperatures below the superconducting transition, as seen in our results. On 
the other hand an energy gap of d-wave symmetry has nodes and thus ensures the availability of scattering states 
even at T = 0. This would produce a quasiparticle damping that drops less rapidly as T decreases below Tc. 

We have also explored the influence of the calculated FL damping on the NMR relaxation rate. We find that the 
damping acts to reduce the magnitude of the Hebel-Slichter peak below Tc, but does not eliminate it. Our results 
are compared to experimental data on three organic superconductors. In the case of (MDT-TTF)2Aul2, our results 
fit the data and reproduce the weak HS peak seen in the experiment. In the case of (TMTSF)2C104, the data show 
no HS peak and deviate from our FL results. We find, however, that if the damping is made anomalously large, our 
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results resemble the data on (TMTSF)2C104, though a small HS peak is still present in the calculations. Since the 
HS peak is sensitive to the BCS coherence factors, the gap symmetry and the quasiparticle damping, the analogous 
"coherence peak" in the microwave conductivity should provide a further test of our analysis. A ppli cations of the FL 
damping to microwave conductivity and surface resistance measurements on the NCCO cuprate |35[ reveal that these 
probes are also sensitive to the self-energy structure |3^ . 

The NMR data on the BEDT organic superconductors exhibit a very rapid decrease below Tc that is similar to the 
behavior seen in the high temperature cuprate superconductors. These cases appear to require additional mechanisms 
to explain the lack of a Hebel-Slichter peak as well as the curvature of the NMR relaxation below Tc- Our microscopic 
analysis yields a relaxation curve that remains considerably higher than the BEDT data in the superconducting state 
even if the Coulomb coupling is artificially increased to unrealistic values. 

The influence of the symmetry of the gap on the relaxation rate may be seen from the behavior of the density of 
states. For an s-wave gap, the superconducting density of states is zero below A and has a square root singularity 
at A. The density of states for a d-wave gap is linear at low energies and has only a logarithmic singularity at A. 
While the Hebel Slichter peak is greatly reduced as a result of replacing a square root singularity by a logarithmic 
one, the relaxation rate in the d-wave case falls off more slowly than the s-wave case as T decreases below Tc- At low 
temperatures, the d-wave relaxation rate varies as , and the BEDT data provide evidence for the latter symmetry. 

We note here that in the case of the NFL damping calculated for a nested Fermi surface by Rieck et al. , the 
quasiparticle damping drops to zero at frequencies lo < 3A more rapidly than the FL damping discussed here. Also 
the nesting generates an enhanced peak in the susceptibility at the nesting vector whose suppression by an isotropic 
energy gap is more pronounced than the case of the averaged susceptibility that dominates the Fermi liquid response 
in the present work. 

The extension of the present analysis to incorporate a d-wave energy gap is in progress. Future studies of a model 
Fermi surface that includes nesting features are warranted in view of the remaining anomalies in the NMR spectra of 
the BEDT organic compound as well as the cuprate superconductors. 
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FIG. 1. The resistivity as a function of temperature is shown for various materials at temperatures above the superconduct- 
ing transition. Note that the resistivity of Pb at these temperatures is linear in T (above the Debye temperature) and is much 
lower than that of the organic metal (TMTSF)2PF6 and the layered compounds TiS2 and Ndi.84Ceo.i6Cu04-j,, whose resis- 
tivities show a quadratic dependence on temperature. This variation resembles Fermi Liquid damping with an anomalous 
enhancement of the electron-electron scattering. 

FIG. 2. The normal state quasiparticle damping r(a;) {i.e. the imaginary part of the self energy) for a two dimensional 
Fermi liquid is plotted as a function of temperature at two different frequencies, ui = 0.01 eV and u> = 10~* eV. The standard 
values of the coupling U N{0) = 1 and bandwidth W = 1 eV produce the small magnitude of F. 

FIG. 3. The momentum averaged BCS susceptibility in the superconducting state is plotted as a function of temperature at 
two frequencies, uj = 0.1 Ao and u) = O.8A0, for a transition temperature Tc = 4.2 K. Note that the overall decrease in {x")q at 
very low T remains a feature at all frequencies, while the coherence peak is only visible at lower frequencies. 

FIG. 4. The calculated quasiparticle damping in the superconducting state is plotted as a function of frequency at throe 
different temperatures. We have used here the parameters UN{0) — 1, W = 1 eV and Tc = 4.2 K. Note that at T = the 
damping vanishes below lo = 3Ao. At higher temperatures, however, there is a finite response below 3A(T) with an unusual 
frequency variation. The damping is quadratic in frequency for uj > 3A(r). 

FIG. 5. The calculated quasiparticle damping in the superconducting state is plotted as a function of temperature at three 
different frequencies. We use the transition temperature Tc = 4.2 K, and UN{0) = 1, = 1 eV. In the normal state, the 
damping is quadratic in temperature. As T decreases below the superconducting transition, the damping drops rapidly to zero 
at low frequencies {i.e to < 3Ao) and saturates at a value that depends on the frequency at higher w. 
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FIG. 6. The NMR relaxation rate 1/TiT data of Takahashi et al. on the organic superconductor (MDT-TTF)2Aul2 are 
plotted as a function of temperature at a frequency w = 42 MHz. The data show a weak Hebel-Slichter peak below Tc = . 

Also shown are the theoretical curves for 1/T\T using the calculated Fermi liquid quasiparticle damping Ffl that includes the 
energy gap (solid line) , and a normal state damping that is quadratic in frequency and temperature at all T and lu (dot-dashed 
line). The dashed line represents the conventional BCS result in the limit F ^ 0. In order to fit the (MDT-TTF)2Aul2 data, 
the amplitude of the Fermi Liquid damping Ffl is increased by a factor of a = 40. The calculated curves as well as the data 
have been normalized to the value of 1/TiT at Tc. 

FIG. 7. The NMR relaxation rate l/TiT data of Takigawa et al. on the organic superconductor (TMTSF)2C104 is shown 

(dots) as a function of temperature at a frequency uj = 42 MHz. This material has a transition temperature Tc = 1.06K and 
does not display a Hebcl-Slichter peak below Tc. Also shown are the theoretical curves for 1/T\T using the calculated Fermi 
liquid quasiparticle damping Ffl where the magnitude of the damping is increased over the calculated value by a factor a. The 
curves are shown for two different values of a. The theoretical curves have been normalized to the value of the 1/T\T data at 
T = Tc. Note that even when a is made unrealistically large, the coherence peak is present in the calculations. 

FIG. 8. The NMR relaxation rate 1/TiT data of Kanoda et al. on the organic superconductor K;-(BEDT-TTF)2CuN(CN)2Br 
is plotted as a function of temperature at a frequency u; = 25 MHz. The solid line represents the theoretical curve for 1/TiT 
using the calculated Fermi liquid quasiparticle damping aFpL where the magnitude of the damping is increased over the bare 
value by a factor a = 40. Also shown is the BCS result (dashed line) for which the coherence peak is more pronounced. The 
theoretical curves as well as the data have been normalized to the value of 1/TiT at Tc = 12 K. The low temperature T^ 
variation of the data resembles d-wave pairing behavior, and the pressure variation of the normal state TiT may indicate spin 
fluctuation contributions which are enhanced by Fermi surface nesting, that is beyond the scope of the present analysis. 

FIG. 9. The total contribution J{E) from the calculated vertex corrections to the NMR relaxation rate is plotted as a function 
of the energy E scaled to the energy gap A at a fixed temperature T = IK. The three curves correspond to superconducting 

transition temperatures of Tc = 100, 20 and 10 K. The correction increases with Tc as soon in the figure. For a fixed Tc the 
vertex contribution at the T = 1 K value chosen here is close to the maximum and falls off as the temperature approaches 
Tc. For Tc = 10 K, the magnitude of the correction at _E = 4A is < 3%, and thus negligible. We also note that since J{E) is 
practically independent of E, the effect of this correction is to renormalize the interaction strength U . 
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